Driving two atoms in an optical cavity into an 
entangled steady state using engineered decay 

Florentin Reiter, Michael J Kastoryano and Anders S S0rensen 

QUANTOP, Danish Quantum Optics Center, Niels Bohr Institute, Blegdamsvej 17, 
DK-2100 Copenhagen 0, Denmark 

E-mail: reiterOnbi . dk 

Abstract. We propose various schemes for dissipative preparation of a maximally 
entangled steady state of two atoms in an optical cavity. Harnessing the natural decay 
processes of cavity photon loss and spontaneous emission, we use an effective operator 
formalism to identify and engineer effective decay processes, which reach an entangled 
steady state of two atoms as the unique fixed point of the dissipative time evolution. 
We investigate various aspects that are crucial for the experimental implementation of 
our schemes in present-day and future cavity quantum electrodynamics systems and 
analytically derive the optimal parameters, the error scaling and the speed of conver- 
gence of our protocols. Our study shows promising performance of our schemes for 
existing cavity experiments and favorable scaling of fidelity and speed with respect to 
the cavity parameters. 
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1. Introduction 

The reliable and efficient preparation of entangled states has been one of the main tasks 
in quantum information science since the birth of the field. The effort has been driven 
on the one hand by the desire to understand these quintessentially non-classical states of 
matter, and on the other by their promise as building blocks for quantum information 
processing tasks. In particular, bipartite maximally entangled states constitute the 
standard of entanglement theory, which in turn is believed to be the main ingredient 
responsible for the additional information processing power of quantum machines over 
classical ones. 

Since the advent of quantum information science, noise has been considered a 
detrimental element in a physical setup, causing decoherence which must at all 
cost be avoided. A few years ago, however, it has been suggested that dissipative 
noise can be used as a resource for quantum information processing, abetting in the 
preparation of entangled states. For instance, Refs. [H El El IH El El E] showed that 
the dissipative dynamics of two atoms coupled to a common reservoir could lead 
to entanglement. These initial ideas were generalized in Refs. [HI El ttH] to show 
that indeed a very general class of states and quantum informations tasks could be 
realized by dissipation. Since then many quantum information processing tasks have 
been reconsidered within the framework of dissipative state engineering. For instance, 
universal quantum computation [8], entanglement distillation and quantum repeaters 
[TT] . quantum memories [12], quantum simulators [121111] and various forms of entangled 
state preparation [ISlIISlIIZlIIHllIHlEOlEIlEaEHlEllE^ have all been 

shown possible using only dissipation as a resource. The physical systems which have 
so-far been proposed or used for dissipative state preparation are: Cavity quantum 
electrodynamics (QED) [Ml [13 Hg EQl |22l [25] , ion traps [21 [131 HI [15] , optical lattices 
P [I9], atomic ensembles [TBI |23l [21], and plasmonic waveguides [211 I2I]. The ffist 
experimental studies along these lines [131 El] have shown these new ideas to be realistic 
and promising as a new path for harnessing the potential of quantum information. 
In this paper, we consider two A-atoms trapped in a single mode cavity QED setup 
[2H1 [2ni [Sni [31] coherently driven by a classical optical field and a microwave or Raman 
field. We demonstrate that a maximally entangled steady state (stationary state) of the 
two atoms can be prepared dissipatively with a very high fidelity by several qualitatively 
different state preparation mechanisms. Which mechanism is desirable ultimately 
depends on the strengths and limitations of the experiment at hand [321 ESI [SI |35l [36] . 
In each of the schemes, the two atoms are rapidly driven into a singlet state, independent 
of the initial state of the system, and without need for any unitary feedback control. 
Consequently, the lifetime of the state is dictated by the lifetime of the experiment. As 
maximally entangled states are an important resource in many quantum information 
processing protocols (ex: repeaters, cryptography), having access to a reliable source 
can not be overestimated. 
Below we will consider a number of schemes: 
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(i) A scheme that employs spontaneous emission mediated by the dark state of the 
atom-cavity interaction. This scheme leads to the highest fidelity of the entangled 
steady state among the presented schemes. In Sec. |4]we analyze this scheme in full 
detail and derive benchmarks for error and speed of the protocol in Sec. [5j 

(ii) In Sec. [6] we present various schemes suitable for existing cavity QED experiments 
which do not have transversal confinement of the atoms, cf. Ref. |32j . 

(iii) For completeness we discuss the scheme of Ref. [22] developed by us previously 
which uses engineered cavity decay (Sec. 11), and an adaptation of a dissipative 
protocol presented in Ref. [20] to A atoms in optical cavities (Sec. [9]). 

For each of the schemes, we identify the relevant interactions by systematically 
truncating the Hilbert space of the problem by an effective operator formalism based 
on second order perturbation theory and adiabatic elimination of the excited states 
[37] . This gives us an effective master equation for the schemes (i-iii), from which all 
of the desired performance measures can be analytically derived. For the scheme (i), 
we analytically derive the optimal steady-state fidelity and the convergence time as a 
function of the system parameters in Sec. |5} In Sec. 10 we analyze its robustness 
against a difference in the coupling of the atoms to the cavity which is present even in 
state-of-the-art optical cavities. In all cases the analytic results are verified by numerical 
simulations. The results obtained are collected and compared in Sec. [7} 
We note that our studies are in a sense related to the more abstract protocols presented in 
Refs. [H El [31 m [5l El [7] , showing that a set of entangled states of two two-level systems 
coupled to a common Markovian bath can be reached by purely dissipative means. 
However, as we aim to demonstrate schemes for concrete cavity QED experiments, we 
focus on maximally entangled states of three-level systems. Our studies can be seen as 
a step towards the desired general tools outlined in Refs. [HI O [10] ■ A highly related 
scheme was also presented in Ref. [25], but a comparison to this scheme is beyond the 
scope of this work. 



2. Cavity QED setups for dissipative preparation of entanglement 

We consider a single-mode cavity QED system consisting of two distantly trapped A- 
type atoms with ground states |0) and |1) and an excited state |e). These levels are 
coupled by a laser, a microwave field or Raman transition, and a cavity mode. Within 
the dipole and rotating wave approximation the couplings of this system are described 
by a Hamiltonian 

H = Ho + Hi,c + -f^laser + -f^MW- (l) 

We assume that the level splittings are the same for the two atoms, and do not 
fiuctuate in time. This can for instance be achieved by cooling the atoms to the ground 
states of identical trapping potentials or by using 'magic wavelength traps'. Then the 
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Figure 1. Cavity QED setup for dissipative preparation of an entangled state between 
two A-type atoms in an optical cavity. Coherent driving U, is performed from ground 
state |0) to the excited state |e). Atom-cavity interaction g takes place between |e) and 
ground state the ground states are coupled by a microwave or Raman transition 
^^MW- Spontaneous emission and cavity photon loss are present as sources of decay. 



Hamiltonian for the atoms and the cavity mode is given by 

H^ = uo,a^a + ^ (u;o|0)j(0| + + a;e|e),-(e|) . (2) 

The couphngs of the levels are expressed by the interaction Hamiltonians 

Aaser = ^6-*^'-* (|e)i(0| + e*<^|e)2(0|) + H.c. (3) 

^Mw = ^ E e-™*|l),(0| + H.c. (4) 
H^c = g^a\e)j{l\+H.c. (5) 

The coherent laser field -ffiaser drives the transition between the ground state |0) and 
the excited state |e) with resonant Rabi frequency VL. The angle determines the phase 
difference of the driving field for the two atoms with respect to the atom-cavity coupling; 
for convenience we assume it on the driving of the second atom. The two ground states 
|0) and |1) are coupled by means of a coherent microwave field or Raman transition Hmmv 
of Rabi frequency I^mw- The atom-cavity interaction H^^ describes that the cavity field 
{a, a)} couples the |1) and |e) states with a strength of g and an identical phase for both 
atoms. Assuming that the ground and excited subspace are coupled perturbatively, the 
system Hamiltonian can be structured according to 

H = Hg + H, + V+ + V., (6) 

with Hg (He) as the Hamiltonian of the ground (excited) subspace and V = V+ + V- 
(V"_ = vl) as the perturbative (de-) excitation term between the ground and excited 
subspaces. We change into a frame rotating at the frequency of the level |0), uq, and 
the frequencies of the laser and the microwave, wiaser and ojmw, to obtain the time- 
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independent couplings illustrated in Fig. [T} 

4 = ^ E + H.c.) + P (7) 

i=l,2 i=l,2 

= A ^ \e)j{e\ + 5a) a + {a)\l)j{e\ + H.c) (8) 

i=l,2 j=l,2 

VV = §(|e)i(0| + e^1e)2(0|), VI = V^. (9) 

Here, A = — Wo~'^iaser and {3 = wi — Wq^i^mw are the detunings of the laser and of the 
microwave field, respectively; a cavity excitation has an energy of 5 = cUc — wiaser + i^mw- 
The noise processes resulting from interactions between the system and the environment 
are assumed Markovian so that the time evolution of the system can be described by a 
master equation p = C[p\ with a Liouvillian in Lindblad form 

p = C[p\ = -t [h, p\+Y. ^^P^l - 2 {^^''^'P + '"^^^'^) • 

k 

The Lindblad operators are associated with the following naturally occurring noise 
processes: (i) loss of a cavity excitation, = a/ko, with a photon decay rate of k, and 
(ii) decay by spontaneous emission from the excited atomic state \e)j into ground states 



|0)j and with Lindblad operators L^aj = v7/2|0)j(e| and L^,i,j = ^/^/2\l)j{e\, 
respectively. Given that the separation of the atoms will typically be more than one 
wavelength for typical experimental conditions we neglect collective components of the 
spontaneous emission. Furthermore, we see from the arguments below that the band- 
width of the laser plays a minor role, as long as it is kept within the linewidth of the 
transition we want to drive, e. g. ~ 6 MHz for the parameters of Ref. |32]. Note that 
for simplicity we assume equal rates 7/2 into the two ground states; an asymmetric 
decay of |e) does not modify the results significantly and can even be used to improve 
the protocols. We ignore the influence from other sources of noise than spontaneous 
emission and cavity decay, such as dephasing due to stray fields or blackbody radiation. 
This is justified if the coherence time of the hyperfine transition exceeds the preparation 
time. As we shall see later the preparation time is on the order of /is justifying this 
approximation for most experiments. 

As a measure for the quality of the cavity we introduce the cooperativity parameter, 

2 

C = —. As we will see later the cooperativity is the main parameter quantifying 
how well the entanglement protocols work. We note here, that we define 7 as the 
decay of population of the excited state and k as the photon loss rate of the cavity, 
which differ from polarization and field decay rates which are also commonly referred 
to as (7, k) [32j by a factor of two each; hence, for the parameters of Ref. [32] we get 
C ~ 17. In this work, we will always assume strong coupling C ^ 1, but we distin- 
guish between the regimes of weak driving {Q, Qmw: P) ^ {.9i 1^1 1) and increased driving 
(fi,fiMw,/3) < (k,7)- 



In the following, it will be convenient to work in the triplet-singlet basis of the 
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Figure 2. Coherent and dissipative interactions between ground and excited states. 
Ground states are coherently excited by an optical field (il) to excited atomic states 
(straight red arrow). Exchange of excitation via the atom-cavity interaction [g) couples 
these to cavity-excited states (wiggled blue arrow). Ground states are coupled by 
a microwave or Raman transition fiiviw- Atomic excitations decay by spontaneous 
emission at a rate of 7 and cavity loss occurs at a rate of n (both indicated with 
dotted arrows). The ground to excited subspace interactions are drawn for (j) — Q 
where excitation happens inside the triplet/singlet subspaces, whereas cj) — n leads to 
crossings between the triplet and singlet sectors. The detunings are defined in the text. 



ground states: {|00), |11), |T), where |00) = |0)i|0)2, |11) = |l)i|l)2, \T) = 

^ (|0)i|l)2 + |l)i|0)2), and \S) = ^ (|0)i|l)2 - |l)i|0)2); the latter is the desired max- 
imally entangled singlet state for all the protocols that we investigate. We further 
define the following excited states which will appear throughout the article: \Tq) = 
^ (|0)i|e)2 + |e)i|0)2), l^o) = (|0)i|e)2 - |e)i|0)2), \T,) = ^ (|l)i|e)2 + |e)i|l)2) and 
l^i) = (|l)i|e)2 — |e)i|l)2). The excited states of the cavity read |00)|1), |11)|1), 
|T)|1), We will truncate the Hilbert space by excluding states with two or more 

excitations (of the cavity or of the atoms), as we always work in the perturbative regime 
with weak driving in our analytical calculations. The states and their couplings (for 
= 0) are shown in Fig. [2] 



3. Effective open system dynamics 



The key to establishing an entangled ground state by dissipative state preparation is to 
identify and engineer decay processes present in the physical system in a systematic and 
reliable way. In the following section we present the methods that are used throughout 
this article to model the open quantum system at hand. We introduce an effective 
operator formalism that allows to reduce the unitary and dissipative dynamics of the 
open-quantum system to the non-decaying ground states and permits us to tailor the 
effective decay processes to achieve a desired steady state. 
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Complex energies and non-Hermitian time evolution of the excited states 

In Fig. [2] we have visualized the coherent and dissipative interactions of the ground and 
excited states present in our system. It is evident that the dynamics of the excited states, 
which are subject to decay, is governed by both, unitary and dissipative time evolution. 
For the excited levels we can combine the real detunings of the levels with imaginary 
terms, that correspond to broadening by decay, to yield complex energy terms. The 
resulting non-Hermitian time evolution of the excited states is expressed compactly in 
terms of a non-Hermitian Hamiltonian 

H^n^H,-lY.^VLk. (11) 

k 

Also referred to as the no-jump Hamiltonian in the language of the quantum jump 
formalism [38], -^nh combines the Hamiltonian of the decaying excited subspace 



with the jump terms of the master equation (10). For the system at hand we find 

i^NH = Ao (|To)(To| + |5o)(5o|) + Ai (|Ti)(Ti| + + 
+ 5o|00)|l)(l|(00|+5i|T)|l)(l|(r| + 
+ 52|ll)|l)(l|(ll|+^i|5)|l)(l|(5| + 

+ (|To)(l|(T| + |5o)(l|(S| + |Ti)(lKll| + H.c.) . (12) 

Here, we have defined complex 'energies' A„ = A„ — ^, 5„ = 5„ — y, A„ = A + m ■ /3, 
and 6n = S + m ■ (3, with m being the number of atoms in state |1) (Aq = A, 60 = 6). 
As will become clear further below, it is useful to set A_i = Ai. 

In doing so, the detunings of the excited states are conveniently combined with their 
decay widths to complex detunings, where both their real and imaginary parts govern 
the strengths of the transitions involving the excited states. 



Effective Hamiltonian and Lindblad operators 

As can be recognized from Fig. |2| the coherent and dissipative couplings of the excited 
states can be concatenated to effective second-order processes between the ground 
states. An example for such an effective unitary process is given by the transition 
|00) —7- |To) —7- |00), facilitated by the coherent driving V of strength Q, resulting in an 
effective shift of ground level |00). In case the coherent de-excitation is replaced by a 
decay an effective dissipative process is formed. Here, |00) ^ |To) ^ |^) is an example 
for effective decay from state |00) to \S), involving spontaneous emission of an atomic 
excitation at a rate of 7. In this manner, all available combinations of weak optical 
excitation, non-Hermitian coupling between the excited states and either subsequent 
weak optical de-excitation or decay bundle together to effective second-order processes 
between the ground states. We assume the optical excitation to be a perturbation 
of the non-Hermitian evolution given of the excited levels by -f^NH- Consequently, their 
population will always be much lower than the population of the ground levels. On 
these grounds, we can perform adiabatic elimination of the excited levels and restrict 
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Figure 3. Propagators in the excited-state subspace. The loop-hke elements A^^g, 
'^n off ^'^^ transition-hke elements g~\f^ contained in H^^^ govern the non-Hcrmitian 
evolution of the excited states. Grouped according to the three interacting and 
two non-interacting excited subspaces these propagators determine the strength of 
effective processes involving the excited states, depending on the state reached by 
initial excitation and the one left by either coherent de-excitation or decay. 



the dynamics to the ground states. 

In Ref. |27] we present an effective operator formahsm based on second-order 
perturbation theory and adiabatic ehmination to reduce the evolution of an open system 
to effective unitary and dissipative processes between its ground states. Applying this 
method simplifies the complexity of the Liouvillian dynamics of Eq. ( 10 ) considerably, 
resulting in an effective master equation in Lindblad form 

p = . [p, ^eff] + 4\p(4'fr)^ - \ {&)^L':,p + p{Ll^YLl^) (13) 

k 

represented by an effective Liouvillian £eff[p] = P- It contains an effective Hamiltonian 
He?t and effective Lindblad operators L^^ 

i^eff = -\V- (i^NH + (^Nh)^) + H, (14) 

4\ ^ L,H^l,V+ (15) 

incorporating the inverse of the non-Hermitian Hamiltonian -f^NH of Eq. (11) and the 
unperturbed ground-state Hamiltonian H^. 

As expected, these effective processes consist of an initial weak optical excitation V+ and 
a final de-excitation or decay depending on their unitary or dissipative character. 
In-between, the inverse non-Hermitian Hamiltonian -ff^H ^^^^ ^ 'propagator', 
representing the non-Hermitian evolution of the excited subspace. We find that its 
elements determine the strength of the effective process depending on which excited 
states take part in it. Its properties will be addressed in more detail in the following 
section. 
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Effective propagators of the excited states 

In the excited-state basis defined earlier, -^nh can be broken up into 5 block diagonal 
elements which evolve independently. -f^NH then also in block diagonal form, and can 
be written out explicitly as 

-^NH = -^NH,|To> + -^NH,|So> + -^NH,|Ti> + -^NH.ISi) + -^NH,|00>|1> 

with three blocks for the interacting excited subspaces 

H-^^^^^^ = A-^g|ro)(To| +5i-lJT)|l)(l|(T| +^-i^(|T)|l)(To| +ff.c.)(17) 
H^i\s,)=K^\So){S,\ + ~5-l^\S)\l){l\{S\+g^^^^ 
H^i\T,)=~^~2lK\Ti){TA + ~5,l^\ll)\^^ 

(19) 

and two blocks for the non-interacting excited states 

^NH,|5i> = ^o,cfrl'S'i)('S'i| (20) 
^NH,|oo>|i> = ^o:efrlOO)|l)(l|(00|. (21) 
In order to keep the notation compact we have defined here 

2 

An,eff = An-1 - ^ (22) 
On 

5n-^ (23) 



■ A„_i 



5'n,eff = Vng - J" , (24) 

are generally complex and their magnitudes have 
ay the role of propagators or complex magnitudes 



The entries of -ff^H' shown in Fig. [s 
dimension of inverse energy. They p 
in the effective ground state to ground state processes of Eqs. (14 15). Each effective 
process that is formed from perturbative optical excitation V+ and subsequent de- 
excitation V_ or decay includes a propagator depending on which excited states are 
involved. 

As can be seen in Fig. sj -ff^n contains both loop-like and transition-like propagators, 
grouped according to the five separable subspaces. Their index n reflects the coupling 
strength between the atomic and cavity excited state of the interacting subspace, the 
latter of which has n atoms in state |1). The states IS*!) and |00)|1) are dark-states of 
the atom-cavity interaction and are uncoupled (n = 0). 

By our definitions of Eqs. (22 23 ) we have associated the loop-hke propagators A^^g and 



'^n off "with the complex detunings of the excited states such that for a vanishing coupling 
g the shifts in A„eff and 5„_eff would disappear, and we would find A^eff = A„_i and 
5n eff = ^n- Similarly, in case of negligible complex detunings A„_i and 5„ the transition- 



like propagators g^^^^ in Eqs. (17 19) would simplify to the inverse of a real coupling 



5'n,eff = \/ng-, the well-known dressed state energy for n atoms resonant with a cavity. 
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All propagators of i^NH and Fig. [s] can be written in terms of a denominator 
Dn = ng"^ — 6n ■ A„_i which equals the reduced determinant of the Hamiltonian of the 
according subspace and is highly dependent on the system parameters. As we will show, 
an appropriate choice of the parameters A, S and /3 engineers certain propagators of the 
excited states, and hence, effective decay processes mediated by them, to be very strong, 
while others are effectively suppressed. Physically this can be understood as shifting one 
of the dressed states into resonance to enhance the coupling. In the following discussion 
we present various schemes that build upon this principle of engineered decay. Here, 
each of the schemes is denoted by the atomic excited state that mediates the engineered 
decay into the desired maximally entangled singlet state \S). 

Applying the discussed methods we will be able to analytically derive the optimal 
parameters and benchmarks for each of the presented schemes, the most important of 
which are the steady-state fidelity with the desired entangled state and the convergence 
time, estimated by the inverse of the spectral gap (as explained later the convergence 
time can be estimated by the spectral gap which is the smallest non-zero eigenvalue of 
the Liouvillian) . We back up all of our analytic results by numerical integration of the 



master equations ( 10 ) and ( 13 ) 



4. A scheme for engineered spontaneous emission mediated by a dark state 

In the following we present a scheme for the preparation of an entangled steady state 
by an engineered spontaneous emission process mediated by the dark state of the atom- 
cavity interaction, l^i). Of the considered schemes it exhibits the lowest error in the 
preparation of the entangled state. 

As for all schemes discussed in this article, we begin by outlining the physical 
mechanisms that underlie the dissipative state preparation, and proceed with a 
discussion of the effective operators. Depending on the driving regime we analytically 
derive the benchmarks for each scheme, in particular steady-state fidelity and speed of 
convergence, from the effective dynamics. 

For the scheme at hand, we choose to engineer the effective decay by spontaneous 
emission into the maximally entangled singlet state to be as strong as possible. To 
this end we set both the optical driving and the cavity transition to be resonant 
(A = 0, (5 = — /3) and the microwave or Raman transition to be slightly detuned 
(/3 7^ 0). Furthermore, we set = vr so that the optical driving crosses the singlet /triplet 
subspace, i. e. |T) ^ l^i) and 15*) 14 |Ti). 

Mechanism of the state preparation 

The mechanism is illustrated in Fig. |4| Population from |T) is excited up to 
\Si) = (|le) — |el)). The atomic excited state IS*!) is the dark state of the atom- 
cavity interaction, ^ac, and is therefore resonant with the optical driving (A = 0). 
Consequently, |T) decays very rapidly over l^i) into \S). On the other hand, population 
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Figure 4. Effective decay processes of tlie sclieme. Tlie optical excitation J7 
drives population from |T) to IS*!). From there it decays spontaneously into the desired 
steady state \S) with a certain probability. As \Si) is the dark state of the atom-cavity 
interaction /fac the initial excitation is not shifted and is close to resonance so that the 
effective decay 7efj from |T) prepares \S) very rapidly. As ifac strongly couples \Ti) 
and |11)|1) with a strength these states form dressed states that are shifted out of 
resonance. Effective spontaneous emission 7eff and cavity loss k^s from 15) into |11), 
mediated by |Ti), is hence effectively suppressed. A microwave/Raman transition (not 
shown) couples |00) and |11) to |T). 



from 15) is excited to |Ti) which is coupled to |T)|1) with a strength of \/2g. This 
strong couphng shifts their dressed states out of resonance by ^\/2g which is much 
more than the natural linewidth. Decay out of \S) is thus strongly suppressed, while 
\T) is rapidly pumped into \S). Accumulation of population in |00) or |11) is prevented 
by the microwave/Raman field that couples the three triplet states |00), |T) and |11) 
and reshuffles population to |T) from which engineered decay prepares 15) again. The 
detuning /3 of the microwave is needed to prevent :^(|00) — |11)) from being a dark 
state of the microwave which would not be reshuffled to \T). The effective processes 
resulting from the coupling are illustrated in Fig. [sja). We note that even though the 
state |Ti) is far out of resonance the desired steady state 15) is still weakly coupled to 
|Ti) by the laser driving; \S) is hence not an ideal dark state. The fidelity of the steady 
state with \S) and the error rate of the protocol depend on the ratio of the rate of the 
dissipative preparation of \S) and the rate of decay from \S). In the following sections 
we will model these processes quantitatively by considering the effective operators to 
derive the optimal parameters and the error of the protocol analytically. 



Effective processes 

We begin our discussion of the effective processes, shown in Fig. [sja), by deriving the 
general effective operators for optical driving with = vr. Given V and the terms 
for the effective processes can be read off directly from the map of propagators in Fig. 
[3| These operators are equally valid for the |5o) scheme in Sec. 11 that also uses (p = n. 
For the effective decay of an atomic excitation |e) into ground state |0) by spontaneous 
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a) |00) b) K) 




Figure 5. Effective ground state processes of the \Si) scheme, (a) Detuning and 
interactions in the shuffling picture, where |00), \T) and |lf) are coupled by i^MW 
to avoid population in |00) and |fl). Engineered spontaneous emission prepares the 
maximally entangled singlet state 15) at a rate of Effective cavity decay out of 

\S) happens at a rate of Kgff. (b) Dressed state picture. Strong spontaneous emission 
oc 7cff reshuffles the dressed triplet states (not shown between \T±)). Population in |5) 
is gained from (oc 7off ) and lost to (oc Kcff ) each of the dressed triplet states. 



emission we obtain the effective Lindblad operators 



-7A{1,2} 
"eff 



± 1 00) (00 1 + v5Z?l (±|T)(T| + \S){T\) + 



+ 



2Ai^eff 

^772^^ 



4Ao,cff 
{\T){S\±\S){S\). 



(25) 



4A2,eff 

The superscript 7 on the Lindblad operators stands for spontaneous emission, for the 
decay into ground state |0), and the index {1,2} for the atom at which the decay occurs 
refers to the upper (lower) set of signs of the terms. Similarly, the effective decay by 
spontaneous emission into ground state |1) is given by 

P7,1,{1,2} _ 



^eff 



(±|r)(oo|-|5)(oo|)±^5^|ii)(T| + 



4Ai,eff 4Ao,eff 
4A2,eff 

The effective decay of a cavity excitation is found to give 



TK 



2fi'2,ofr 



fi'l.cff 



\s)m. 



(26) 



(27) 



Finally, the effective unitary processes are given by 



eff 




1 00) (00 1 - Re 



\S){S\+H,. 



4A, 



O.eff 



\T){T\ 



(28) 



where Re[ ] denotes the real part of the argument. While the above effective operators 
hold whenever = tt, we can simplify them for the particular scheme at hand by 
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discussing the propagators for the parameter choices made in the previous section. In 
the absence of an atomic detuning, A = 0, the complex energy of l^i), as the dark- 
state of the cavity interaction, is given by Ai = /3 — y. As we will discuss below it is 
desirable to have /3 oc f2 so that for the assumption of weak driving (f2 ^ 7) we can 
write Ai ^ — ^. The propagator of the effective IS*!) -mediated decay processes then 
simplifies to 

(5i|^NA|5i) = Ao,eV^-|-. (29) 

Hence, the effective decay processes mediated by l^i) that incorporate this propagator 
are tailored to be very strong compared to the decay out of the singlet state which 
involves the subspace consisting of the states |Ti) and |11)|1), and the transition-like 
propagator 

(l|(T|^^Al^l)=^2":eff^^- (30) 

The last denominator reflects the strong shift of the dressed states of |Ti) and |11)|1) 
out of resonance, slowing down the effective decay out of 15"). Consequently, we have 
reached IAq ^jjI ^ |^„^,effl that effective processes mediated by the dark state l^i) are 
engineered to be much stronger than those involving other states, in particular |Ti). We 
have thus found that the triplet ground state |T) undergoes rapid effective spontaneous 
emission at a rate oc I/7, while cavity decay from \S) oc k/q"^ is suppressed in the strong 
coupling regime where C ^ 1. For finding the steady state we can drop the suppressed 
terms unless they affect the singlet state. In addition, the spontaneous emission from \S) 
oc '-yn? / g'^ is negligible and will be ignored. The effective decay processes then simplify 
to 

Ll§^^^^^^ = ±^^\T){T\+^^\S){T\ (31) 
L^^'^^-^> = ±zv^|ll)(T| (32) 
4\ =^\U){S\. (33) 

Here, we have set = |(11|4\I'5)P = ^ and 7eff = \{S\Lls''^'^'^\T)\^ = f^. 
Furthermore, for the scheme at hand the effective Hamiltonian ifgfr of Eq. (28) is 
well approximated by the unperturbed ground-state Hamiltonian Hg 

H, = ^ (|00)(T| + |T)(11| + H.c.) + /3 (2|11)(11| + \T){T\ + \S){S\) , 

where we have neglected the minor effective shifts (9(fi^). The resulting effective decay 
processes of this scheme are illustrated in Fig. |5]a) together with the microwave /Raman 
reshuffling. The singlet state IS") is efficiently prepared from |T) by spontaneous emission 
at a rate of 27eff. The singlet \S) decays by effective cavity loss k^s into |11). The 
mechanism that allows us to engineer a strong effective spontaneous emission from |T) 
into IS*) at the same time causes strong decay at a rate of -ij^s from |T) into |11). Hence, 
accumulation in |11) needs to be inhibited by the microwave/Raman shuffling ^mw- The 
fidelity of the steady state with the desired entangled state will be be derived analytically 
in the following sections after changing into a dressed state picture with respect to fluw- 
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Parameter analysis at weak driving 

We first analyze the dynamics of tliis l^i) sclieme for weak optical driving VL. After a 
basis transform into a dressed ground state picture, this assumption allows us to reduce 
the dynamics to rate equations for the ground state populations. From these, we derive 
the important benchmarks for a comparison of the presented schemes; the steady-state 
fidelity with the desired entangled state, and the spectral gap as a measure for the rate 
of convergence. 

The basis used so far, involving the triplet states |00), |T) and |11) coupled by ^^mW) 
will be referred to as 'shuffling picture' in the remainder of the paper. We now simplify 
the analysis by identifying a basis in which non-diagonal elements of the density matrix 
of the reduced system are suppressed, as a consequence of the weak driving. It is then 
possible to express the dynamics as a set of linear rate equations. The basis of the 
new 'dressed state picture' contains the original singlet state \S), and the three dressed 
triplet states 

|r±) = -1/2(5 T 1)|00) + 1/2(5 ± 1)|11) + A/V2\T) (34) 
\Tr) = A/^/2\m)-A/V2\ll) + B\T), (35) 

where we have defined A = ^mw/ ^/^mw + /^^ ^^'^ B = (3/ -^/fi^^ + (5"^. In this basis 
the effective Hamiltonian is diagonal 

^eff = 5^ (/3 ± (5/3 + A^Mw)) |T±)(T±| + (3 {\Tr){Tr\ + \S){S\) . (36) 

+/- 

The parameters A and 5, and hence the ratio of /3 to f^Mw, determine the contribution 
of |T) to each of the dressed states. We find that the optimal fidelity is obtained at 

A = ^^J^, B = and /3 = f^Mw/v^- Here, each of the dressed states contains an 
equal share of the triplet state |T) from which \S) is prepared. 

In the weak driving regime (fi ^ 7), the rephasing of the dressed states is much faster 
than the effective decay 7efr oc — . Consequently, in the new basis the evolution of the 
coherences can be dropped from the master equation. The dissipative time evolution is 
then well approximated by a set of coupled linear differential rate equations. The rate 
equation for the population of the singlet state P5 is 

Ps = ^ {Pt, + Pt^ + Pt,.) - -^Ps. (37) 

Here, we have used that for weak driving the decay from the three dressed triplet states 
into the singlet are of the same strength due to the equal weight of |T) in the dressed 
states and equals one third of the total rate Q/i'j. 

Derivation of the static error 



From the rate equation (37) we can derive the fidelity of the steady state with respect 



to the maximally entangled singlet state as Fs = lim Ps- Equivalently, the error of the 
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protocol is found as the stationary population of the undesired triplet states (1 — 



s 



Inserting 1 — P5 = Pj^^ + P^^ + i^ito Eq. (37), we use = and obtain for the 
static error of the protocol 

37K _ 3 



with the cooperativity C as defined in Sec. [2j Eq. (38) indicates that in the strong 
coupling regime (C ^ 1) the only non- negligible error term scales linear in C~^. This 
linear scaling of the error in the cooperativity is similar to Ref. [22], but the constant pre- 
factor is improved from | to |, which could be important for experimental realizations, 
as discussed in the comparison section [7| 

Derivation of the spectral gap for weak- driving 

The quality of a dissipative state preparation protocol is determined by two main 
benchmarks: (i) fidelity of the stationary state, and (ii) speed of convergence of the 
protocol. We now consider the latter. Estimating the speed of convergence is in general 
a difficult task, but, for small systems, the spectral gap of the LiouvilliaE[|] is a very 
good estimate of the rate of convergence. The spectral gap can be seen as the decay 
rate of the slowest-decaying quasi-stationary eigenstate. If the gap is small, then the 
quasi-stationary eigenstate remains populated for a long period of time, whereas if the 
gap is large, then all eigenstates except the stationary one get depopulated rapidly. 
In the setting at hand, the gap can in fact be read off from the expressions for the 
effective decay process L^g'^^'^^ which have a rate 7efr = As stated above, the 
dressed states each contain an equal share of :^|^)- Hence, the singlet is prepared 
equally fast by decay of the three dressed states at an individual rate of 2 ■ | ■ which 
results in the spectral gap 

' - h <^«' 

Furthermore, A is recognized as the eigenvalue for the according lowest-lying eigenvector 
of the Liouvillian Pj- = | (Pt+ + -Pr„ + -Pr^), as can be seen from Eq. (37) above. 

To confirm these predictions we have performed numerical simulation of the dynam- 
ics of the system. Fig. M (a-b) summarizes the results obtained by numerical integration 



of the master equation of the full system in Eq. (10), consisting of ground states and 



singly excited states. These curves are plotted together either with those from the effec- 



tive master equation in Eq. (13) in (a), or from the rate equations (such as Eq. (37) for 
the singlet) in (b): For a weak driving (fi = :^) we compare the numerically obtained 
curves of the population dynamics of the full master equation (solid lines) with (a) 
the effective master equation in the shuffiing picture and (b) the rate equations in the 

i The spectral gap of a Liouvillian £ is defined as the magnitude of the smallest (in absolute value) 
non-zero real part of the eigenvalues of C, where the Liouvillian is written as a linear operator in the 
matrix units basis, see Ref. [32| for further details. 
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Time (1/g) Time (1/g) 

Figure 6. Evolution of tlie system towards the entangled steady state for (7, k) = 
if) similar to Ref. corresponding to C w 17. Dynamics of the full master 

equation (solid lines) are compared with effective dynamics in the shuffling picture (a, 
c) and rate equations in the dressed state picture (b, d) (dashed) for weak-driving 
(a-b) and increased driving (c-d). Starting from a completely mixed triplet state (see 
legend for details) the system evolves towards the maximally entangled singlet state 
(jS) - green) approaching the steady-state fidelity (dotted line, indicated). For a weak 
driving oi — {'•^~^) the dynamics are completely described by rate equations of 
the populations, shown in (b). At increased driving = ^ (c-d) the dressed effective 
operators (long dash in c) are almost indistinguishable from the full dynamics, while 
simple effective operators and rate equations (short green dash in c), d) resp.) exhibit 
increasing inaccuracies. For all curves the optimized parameters 51 — 2^/''f2MW — 2^/'*/? 
(see also Sec. [s]) were used. 

dressed state picture (both dashed). We see that in this regime the full and the effective 
dynamics of the system are in excellent agreement. In addition, the analytical quantities 
derived from the rate equations in this Section are found to describe the fidelity of the 
steady state with the maximally entangled singlet state and the convergence time very 
accurately. 
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5. How fast can two atoms be entangled by dissipation? 

The fidelity of the prepared state with respect to the desired state gives us a measure 
of the quality of our scheme once the system has reached equilibrium. The second 
figure of merit of a dissipative state preparation protocol is the time required to reach 
convergence. In this section, we analyze how fast the scheme presented above can be 
performed. We emphasize in particular the trade-off between the speed of the protocol 
and its fidelity. 

Speeding up the state preparation can be done by increasing the optical driving. Indeed, 
as can be seen from Fig. |6] (c, d), using an increased optical driving Vt improves the 
convergence time by several orders of magnitude at the expense of only a few per cent 
additional error. The main reason for the decrease in fidelity is that the strong driving 
requires a strong microwave shuffling of the population of the triplet ground states. The 
microwave field, in turn, shifts the ground states out of resonance. This results in a 
decrease of the fidelity at increased optical driving which we will investigate in detail 
below. 

In order to model the effective dynamics of our scheme accurately even for increased 
optical driving, we begin this section by introducing an extended effective operator 
formalism to account for the coherent coupling of the ground states which has so 
far been ignored when adiabatically eliminating the excited states. We then proceed 
to analytically derive the scaling of the two main performance measures, error and 
convergence time (spectral gap), as a function of the strength of the coherent driving, 
and perform a study of the optimal preparation time of an entangled state of a given 
fidelity. 



Effective processes in the presence of ground state dressing 

In the following section we will discuss an extension of the effective operator formalism 
presented in Section [3j It allows us to include the effects of increased driving and 
dressing of the ground states and to derive the dynamic benchmarks of the scheme at 
hand. 

So far, we have worked within the weak driving limit {Vt ^ 7, k), where simple 
perturbation theory holds very reliably. We now want to consider how our scheme 
behaves as we approach the increased driving regime. From the curves of Fig. [6] a) 
and b) we notice an excellent agreement between the dynamics simulated with the full 
and effective master equation, and the rate equations for a weak optical driving as low 
as VL = On the other hand, for Fig. |6] c) and d), we have used Vt = ^ which is 
clearly beyond the weak driving limit. Here, the previously employed simple effective 
operators and rate equations become increasingly inaccurate. This is due to the fact 



that in our derivation of Eqs. (14 15) we have neglected the influence of the ground- 
state Hamiltonian on the effective processes. As we derive below, Vtuw has to be 
proportional to Vt so that an increased f2 also leads to a higher I^mw- For the case 
at hand, this influence induces a shift of the ground states by the microwave driving 
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('dressing'). In a frame where is diagonal we can include these effects by applying 
more general effective operators 



Hp(r = 

2 



H, (40) 



4\ = ^E(^NH-i5/)"'v:f, (41) 



where Ei is the energy of the initial ground state / and the excitation from it. For 
non-negligible ground state energies Ei, the elements of are generally replaced by 
-^NH ~^ {Hnh — Ei)~^. Yet, we note that in order to capture the effects of the dressed 
ground states it will not be necessary to keep the dressed state energies Ei in all propa- 
gators, but only those engineered to be strong. In fact, numerical curves obtained from 
these extended operators match the evolution of the full master equation very accu- 
rately, as can be seen from Fig. |6]c). 

As soon as ground state dressing is taken into account, the decay rates from the dressed 
triplet states \T±) and \Tr) into 15*) are no longer identical. Effective spontaneous 
emission from the dressed triplet states into the singlet mediated by l^i) exhibits the 
non-degenerate propagators 



{Si\ (^Hnr — Et^ \Si) = ^Ao,eff — Et± 
{Si\ (JInu — Et,^ \Si) = (^Ao,cfr — Et^ 



-1 




resulting in detuned decay rates '"f^g^^ = 7eff ■ 7^/(7^ + ^mw)) "while 'j'^g' = 'jes and 

(S) 

nlff ~ Kes are effectively unchanged. Taking these into account we can set up the rate 
equations in the same manner as in the weak driving case of Sec. |4j 

Derivation of the error and of the spectral gap in the presence of ground state dressing 

Employing the state-dependent decay rates 7^^ ^ the additional error originating from 
the dressing of the triplet states is derived the following way: 

Despite the different decay rates into the singlet, the population of the dressed triplet 
states is kept close to an equilibrium by strong dissipative shuffling oc I/7 in-between 
them (see Fig. [s] b). Consequently, an equal mixture of the triplet states Pj- = 
I {E't+ + -Pr_ + Pxr) is the slowest decaying eigenvector of the Liouvillian. Using this 
definition we set up the rate equation for the population of the singlet state 

Pq = — ^ Pt P^ f44) 

127 7^ + 6^^"^ 8^2^^- 

While the loss of population from the singlet by cavity decay Keg is unaffected by the 
dressing, the decay rate of the triplet population through spontaneous emission has 
now become dependent on f^MW^ Introducing Pq- results in an effective decay rate of 
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Figure 7. Fidelity and spectral gap as a function of driving strength, (a) Fidelity 
of the steady state with the singlet state and (b) spectral gap as a measure of the 
speed of convergence with respect to the optical driving strength £7. Analytical results 
(blue, short dash) are in very good agreement with numerical curves obtained from the 
full (dark green) and effective (green dash) Liouvillian even at « 7/2 « 0.2(7. For 
all curves the cavity parameters (7, k) = (^, |f) [32] j corresponding to C « 17, and 
the optimized driving parameters Vt — 2^/'*nMW = 2^/^/3 were used. Note that the 



analytical curve in a) contains terms that are not included in Eq. ( 46 ) (see discussion 



ofEq. (56)) 



7efr 



\ ' V ' 7cfr with a factor 77 



7gg'*. For the steady state (P5 = 0, Ps 

3 72 + 6fi 
2C ' 



2 

MW 



2 + 902 



-^•5) Stat 



1 + 



402 

7^ 



(45) 



originating from averaging the decay rates 
1) we derive the error 
3 

1 - ^5)d.es • ' (46) 

As can be seen from the second step where we have expanded for Omw ^ 7, the errors 
decouple into the static error, derived in Sec. |4| and another dynamic error (1 — Fg)^^^^. 
The latter emerges from the dressing of the ground states by Omw- Just as the static 
error, it decreases linearly with one over the cooperativity C~^. 

In the same manner, the spectral gap in the presence of ground state dressing is found 



from Eq. (44), determined by the decay rate of the lowest lying eigenvector P7-, 

02 



A 



7 



2 + 202 



127 7^ + 60^^- ^^^^ 
This result can also be derived more rigorously if we set up the full rate equations and 
extract the spectral gap as their smallest non-zero eigenvalue. We find 

02 (572 + 1802,^ - v/97^ + 847202,^ + 3240t 



A 



"MW 



247(72 + 602,^) 



„2 I 002 



"2 + fi02 



(48) 
(49) 



_ o^ r 
~ 12^ ' ^ 

In the last line we have used 7 ^ f^MW and expanded up to second order in Omw which 
reproduces the result of Eq. (47). For Omw the derived expressions reduce to the 
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weak driving case as expected. In Fig. [7]b) we have plotted the analytic result for the 
spectral gap (from Eq. (48)) with respect to the optical driving fi, together with the 
numerically obtained spectral gap of the full and effective Liouvillians of Eqs. ( 10 ) and 
(13). We see that also for increased driving the curves are in good accordance. 



Beyond rate equations 



So far, we have carried out our analytic study using rate equations formulated in a 
dressed state picture, where is diagonal. In the same picture, we have included the 
dressed ground state energies into the effective operators. However, in order to fully 
describe the system, in particular the effects at low microwave driving, we change back 
into the original 'shuffling picture' with triplet states |00), |11) and |r) coupled by 
f^MW / v^- Introducing new decay rates we can write the dressed effective operators as 



^7,0,{l,2} 
^eff 



f7,l,{l,2} 

JK 



±t^,\T){T\+i^,\S){T\T 

TXa\T)m-Xa\S)mTxl\T){n\ 

Tv^Xa|ll)(00|Tv^x:|ll)(ll|±M 
+ v^|ll)(5| -2y^|5)(00| 



-Z\s){n\ 

^d\ll){T\ 



(50) 
(51) 
(52) 
(53) 



with the previous but shifted effective spontaneous emission rate 7^ 
7eff • r] and an additional spontaneous emission process activated by f^MW, with an 



87 72+6CI 



amplitude Xa 



Here, Xa carries a phase; the according effective decay 



2^ 7^+6% 

rate is defined by 7a = IXaP, with 7a ^ 7d- The effective cavity decay Kefr = |p- as the 
main loss process remains unaffected by the dressing since VLm\n ^ 9- The shifts of the 
effective Hamiltonian are again negligible so that i^efr ~ -^g- 




b) 
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Figure 8. Effective processes at increased driving, (a) Unaffected (black) and reduced 
(blue) processes. A detuning /3 = I^mw/v^ and dephasing 27^ of state \T) retard the 
recycling of population |11) \T) — > \S). (b) Effective decay processes activated by 
dressing of the ground states by I^mw (red) . 



The effects of increased driving are illustrated in Fig. [8] The dressing of the triplet 
ground states \T±) by ^mw causes a reduction of the advantageous decay processes to 
7rf from 7eff, while the amplitude of the detrimental process Keg remains unchanged. 
The scaling of 7d with rj is the result of averaging over the decay rates of the dressed 
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triplet states through back-transform. In addition, new decay channels at rates of 
7a are activated for high J^mw and increase the error of the protocol by accumulation 
of population in state |11). 



Derivation of the recycling error and optimal reshuffling 

Above we have derived the static error and the driving-dependent error originating from 
the shift of the ground states. An additional error emerges from the dynamics of the 
coherences which is not caught by rate equations of the populations: 
From Fig. |6]a), c) we can see that accumulation of population in |11) is the bottleneck of 
the scheme. Coherent reshuffling fijviw from |11) and |T) is used to recycle the population 
of |11). Hence, the additional error of accumulation of |11) is more pronounced the 
weaker Omw is compared to Vt, regardless of the driving regime. The effective dephasing 
of |T) at a rate of 7d and the detuning of /3 = f^Mw/v^, however, tend to slow down 
the recycling process. Still, justified by its rapid decay of 67^ altogether, the state |T) 
can be considered transient and can hence be adiabatically eliminated. In addition, we 
can ignore the evolution of |00) and the much weaker effective decay processes activated 
by l^Mw which have rates 7a. After adiabatic elimination of the rapidly dephasing 
coherences (pii,r, Pt,ii ~ 0) and the population of |T) (pt,t ^ Pii,ii) we can write the 
dynamics in terms of two rates affecting the population of the desired singlet state 

P,^-«:effP5+^^%Pii. (54) 

Thus, for the steady state (P5 = 0, ~ 1) we derive the error 

(1 - Fs),ecy - -F^2 • (55) 



In order to make sure that the errors of Eqs. (45) and (55) are actually sufficient to 
describe the fidelity of the protocol at increased driving we also derive the steady state 
from the Liouvillian dynamics. To this end we solve the master equation represented 
by the effective Liouvillian Ces for pj^k = for all j, k. 

The extended decay rates 7d and 7a also hold for stronger driving. Since the shuffiing 
fluw is still much lower than the spontaneous emission 7, we also neglect dephasing 
originating from the f2MW-activated processes (7a), as well as dephasing at a rate of Kes, 
in the presence of dephasing at rates of 7d (7a5 ^^eff ^ 7d)- Normalizing the obtained 
expression for the steady state and expanding it up to the second order in Q and fluw, 
we extract the complete driving-dependent error as 

" "^'^-^ + ^''^ 

= (1 - ^S).tat + (1 - ^5)drcs + (1 - Fs^^cy ■ (57) 



This is exactly the sum of the driving-dependent errors of Eqs. (45) and (55), expanded 
for small ^mw- We see that in fact the errors decouple. As one of these terms scales as 
fi^^ and the other as ^mw^ optimum for ^mw is a trade-off between fast recycling 
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requiring large I^mw and the need not to shift states out of resonance favoring small 



Omw- We use the result of Eq. (56) to derive an optimal reshuffling of 



^Mw.opt - (58) 



Inserting r2MW,opt into Eq. (56) we obtain the combined error 

(l-^.)comb=|.(^l + v^(^)')- (59) 

We will use this result below to discuss the scaling of the error with the speed of 
convergence and with the preparation time. 

In Fig. [T] a) we plot the analytical result for (1 — -Fs')comb together with numerical 
curves obtained by extracting the steady state from the full and the effective Liouvillian 
for different optical driving fi, using the parameters of Ref. [32], (7, k) = (^5 §2)5 
corresponding to a cooperativity of C ~ 17. Note that for the analytical curve in Fig. 



7|a) we have not discarded terms of higher order in C~^, as in Eq. (56), but have kept 
terms up to second order in after solving for the steady state of the Liouvillian. 
For higher cooperativies, the higher orders become negligible and the expression for 



(1 - Fs)^„^^ reduces to Eq. (56). 



Fig. [7]b) contains the analytical and numerical results for the spectral gap, as a measure 
for the convergence rate. We find that for our initial assumption of weak driving [Vt <^ 7) 
the analytic results for the scaling of both important performance measures, error and 
spectral gap, with the driving strength are very accurate. In addition, we find very good 
agreement with numerical results obtained from both the full and the effective master 
equation even up to an increased driving of Q ^ 7/2 ^ 0.2g. Despite the increased 
driving, the population of the excited states, in particular the close-to-resonant IS*!), 
does not exceed ~ 5% for i7 = 7/2 (cavity parameters of Ref. jS^) so that both the 
initial truncation of the Hilbert space to ground states and singly excited states as well 
as the concept of the effective dynamics of the ground states are well-justified even in 
the regime of increased driving. 

Performance of the scheme at increased driving 

We evaluate the performance of the scheme, this time at increased driving, by estimating 
the trade-off between fidelity and convergence time. To this end, we use the results for 



the driving- dependent spectral gap and error of Eqs. (48) and (59) from which we 
eliminate the driving Q. In doing so we find 

3 2^27 + 21A A<7 3 18^/2KX 
2C2y27 - 27A 2C ^ ~^ 
For strong coupling g ^ (7, /«), or sufflciently high cooperativities C ^ 10, the static 
and the dynamic error decouple when the expression is expanded in A. Thereby we 
obtain the slope of the tangent of Fs for a small spectral gap A. The analytic result 
shown in Fig. M agrees very well with the numerical results obtained from the full 
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and dressed effective master equation as long as the assumption of perturbative optical 
driving is justified. For very rapid state preparation, the analytic expressions reproduce 
the decreasing trend of the numerical curves correctly. In Sec. [7] we compare these 
benchmarks for the scheme at hand with the ones for the schemes presented in Sec. [6} 
Sec. lUand Sec. [TTl 




Figure 9. Performance of the dissipative state preparation at increased driving, a) 
Tire fidelity of the steady state is lowered by the increase of the dynamic error when the 
spectral gap increases, (b) Error of the protocol (left axis) and optimal driving strength 
(right axis) vs. desired preparation time. Analytic results (blue, short dash) are in 
good agreement with numerical curves obtained from the full (dark green) and effective 
Liouvillian (green dash). The cavity parameters (7, k) = (^iff) [32] (C « 17) and 
the optimized driving parameters = 2^/'*i7MW — "27^^(3 were used; in (b) we also use 
and plot the optimized optical driving ilopt (red dot) of Eq. (62 1. 



Scaling of the dynamic error with the preparation time 

In the discussed setting, the scaling of the error and spectral gap provides an estimate 
of how fast the population decays into a desired steady state and to which extent the 
fidelity is lowered by an increased driving. For preparation within a fixed time one will 
thus have to make a compromise between convergence rate and the detrimental effects 
of increased driving. These two effects can be used to derive the optimal driving for 
a desired preparation time. To this end, the error of the protocol with respect to the 
preparation time t, consisting of a static and a dynamic part, can be written as 

{l-Fs){n,t) = ^ + fn' + le-''''/\ (61) 

with / and r specified below. Here we have assumed that the evolution begins in a 
complete statistical mixture of the four states. We minimize the error for a given state 
preparation time t by taking its derivative with respect to and obtain for the optimal 
driving strength 

^opt = -7log^- (62) 



Driving two atoms in a cavity into an entangled steady state using engineered decay 24 



Above we have derived the combined static and dynamic error and the spectral gap (Eqs. 



(59) and (48)) which we now associate with the quantities / and r. With / 



3k 



V2r 



7 




and r = I27 we get fr = We thus obtain for the optimized error of the protocol 

• (63) 

In the second line we have expressed the optimized error in terms of the cooperativity 
and the ratio of the decay rates 7 and k. We find that - apart from the linear static error 
(^) ~ the above expression for the optimized error also exhibits a favorable scaling of 
the dynamic error part with the square-root of the cooperativity. We plot this analytical 
result in Fig. |9]b), together with curves obtained numerically from the full/effective 
Liouvillian, using the cavity parameters (7, k) = ||) [32] (C ~ 17) and the optimized 



8 ' 32^ 

driving parameters = 2^/'^f2MW = 2^/''/3, as well as the optimized optical driving fiopt 



of Eq. (62). We find good agreement even for reasonably short preparation times 



lO^g ^ where we get fidelities above 90%. 



6. Schemes for various experimental situations 

The 15*1) scheme for effective spontaneous emission mediated by a dark state, discussed 



in the preceding sections, as well as the \So) scheme, discussed in Sec. 11, both assume 
a phase difference of = vr between the optical driving of the two atoms. While in 
present-day cavity experiments the position of the atoms along the cavity axis with 
respect to the cavity standing wave is well-controllable within the Lamb-Dicke regime, 
their transversal position is not necessarily confined. If the atoms are driven by laser 
fields oriented transverse to the cavity axis this results in a random phase factor e*'^''''-*'* 
(with wave vector k and relative position of the atoms r{t)). Hence, the assumption of a 
relative and stable phase relation (p rules out common transversal driving of the atoms 
in the absence of transversal trapping. 

This section deals with better suited alternatives for today's cavity experiments: We 
present a |To) and a |Ti) scheme that can be implemented with a driving with = 0. 



This phase relation can be obtained by driving the cavity with a strong laser which 
is strongly detuned from a cavity mode but near resonance with the atomic transition 
|0) -H- \e). The detuned drive creates a coupling mediated by the off-resonant cavity 
mode and the phase relation will then be set by the phase of the cavity mode. If the 
cavity driving mode and the mode used to create the entanglement are commensurate, 
this will ensure that we have the phase relation = 0. 

As an alternative, we also discuss the possibility to use common addressing by a trans- 
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verse laser at an arbitrary relative phase 0, which results in an effective combination of 
the \Tq) and l^o) scheme. 

A scheme for engineered decay mediated by \Tq) 

For the scheme at hand, we choose to use the subspace containing the atomic excited 
state |To) and the cavity-excited state |T)|1), in order to realize strong spontaneous 
emission from |00) into \S), mediated by |To). For this |To) scheme we use non- 
vanishing laser (A) and cavity (5) detuning, but a vanishing detuning /3 = of the 
microwave/Raman field. This means that A = A„ and 6 = 6n for all n. 



\s,)=^{M-\ei)) |r„)=J=(|o.)+|.o)) 



r 1 



n 



1^) 



r r 



n 



00) 



b) 




1^) 



Figure 10. Mechanism and effective processes of tlie |To) scheme, (a) Dissipative 
preparation of the maximally entangled singlet state 15*). For an appropriate choice of 
atomic and cavity detuning A and S the cavity interaction of strength V^g enhances 
the effective spontaneous emission jcB from state |00) by shifting the lower dressed 
state of \To) and into resonance with the optical driving H.. States |T) and 

1 11) are coupled to |00) by a microwave field or Raman transition riMW (not shown). 
Effective decay from \S) is suppressed as it involves the atom-cavity interaction dark 
state IS*!), the detuning A of which is not compensated, (b) Effective level scheme and 
ground state to ground state processes for the presented scheme. 



Mechanism of the state preparation The working principle is illustrated in Fig. 10 
Population from state |00) is excited to \Tq) by a weak optical field of strength fl. The 
atomic excited state |To) is coupled to the cavity excited state |T)|1) by the atom-cavity 
interaction f^ac- Due to the strong coupling (g), the states |To) and |T)|1), initially 
detuned by A and 6, form dressed states. Treating the detunings of the excited states 
as complex, as discussed in Sec. [3} the energies of these dressed states can be written 
as 



E^ = ^±l±U(A + 6y -4{A6~g^], (64) 



2 2 

where A = A — y and 5 = 5 — y . We engineer an efficient spontaneous emission process 
that prepares the singlet state \S) by setting the cavity detuning equal to the cavity line 
shift S = g'^/A. With this choice the lower dressed state of \Tq) and is shifted 
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close to resonance, Re(i?_) ~ 0. Consequently, population from |00) - which is coupled 
to \Tq) - is rapidly transferred to IS*) by spontaneous emission via the lower dressed 
state of |To) and |r)|l). On the other hand, decay out of \S) involves excitation of IS*!). 
Since 15*1) is the only dark state of the atom-cavity interaction, its detuning A means 
that it is not in resonance. Hence, the decay into the singlet state IS") is engineered 
to be much stronger than the decay out of \S) so that the maximally entangled state 
15) is efficiently prepared. The atomic detuning provides a trade-off between virtual 
character of the excited dressed states on the one hand, and spontaneous and cavity 
decay on the other hand; by setting A = g^/^ we minimize the decay width of the 
dressed states, E_ ^ lm(£'_) ^ |(Ak + (57). Furthermore, coherent coupling of the 
triplet states |00), |11), and |T) by the microwave/Raman field ^mw guarantees that 
all triplet states decay rapidly towards the singlet state \S). 



Effective processes For the discussed scheme we have assumed the relative phase of 
the optical driving between the atoms to be zero (0 = 0). Using this, the effective 
processes can be read off from Figs. [2] and [3j Given our choice of the parameters 5, A 
and /3 the terms of the effective detunings can be simplified to obtain the scaling of the 
effective decay processes. Using ^i^gg i^/fn, ^2,eff ~ 5'/v^, Ao,cff ~ A, Ai^cfr ~ -^7 
and A2,cff ~ —A the effective operators simplify to 

Ll§^^'''^ = ^V^mm ± V^\T){S\ (65) 

L'J'^''"^ = ^v^mm T \s)m) ± v^\n){s\ m 

4\ =-i^\T){00\. (67) 

Here, the spontaneous emission processes L^g'^^'"^^ transfer population from |00) into the 
desired state jS") at a strongly enhanced rate of 27+ = Loss from the singlet state 



also occurs by spontaneous emission at a much weaker rate of 7_ = As opposed to 



2 



the \So) and l^i) schemes, the effective cavity decay, here with a rate of = f^, does 



not directly affect the singlet state. The effective processes are illustrated in Fig. 10 b). 



Parameter and performance analysis Setting up the rate equations in the same manner 
as for the previous scheme is straightforward. We obtain for the error and spectral gap 

(l-^^)|To)=^ (68) 

A|To> = (69) 

Both the error and the spectral gap are found to have the same scaling with the 
parameters of the system as the \Si) scheme, but exhibit different proportionality factors. 
The performance of the schemes presented in this section is optimal for an Quw in the 
interval of ^mw = f to I^mw = f ; the latter value is used for the simulations below. 
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Laser driving with random relative phase between the atoms: the |To)/|5'o) scheme 

One can conceive of experimental situations for which neither transversal confinement 
of the atoms (and hence a stable phase relation), nor cavity driving are available. In 
the following we argue that using laser addressing of the two atoms at random relative 
phase can be suitable for the preparation of an entangled steady state of high fidelity. 
Apart from the driving phase the conditions for the operation of the \Tq) scheme 
of the previous section and the 15*0) scheme presented in Ref. [22] which is briefly 
discussed in Sec. 11 are identical; in particular 5 = g'^/A. Transversal drift of the 
atoms results in a random with fluctuations much slower than the couplings of the 
system. Depending on the actual value of the driving either crosses the singlet /triplet 



subspace, |00) \Sq), as illustrated in Fig. [l4|or stays within the subspace, |00) ^ \Tq) 



n 



similar to Fig. [TOj Therefore, effective decay channels are instantaneously weighted with 
0, as ■jesit) = (7cfr)|ro) " cos^ + (7cff)|5o> ' ^^^'^ Thus, the system mechanisms are 
an combination of the two individual schemes. Averaging the decay rates results in a 
combined error and spectral gap 

(l-^^)|To)/|5o) =^ (70) 

(9 - 2v/3 - v/5)^]^ 

A|To)/|5o) = ^ (71) 

We conclude that a setup with arbitrary driving phase is also suitable for an experimental 
realization of a high-fidelity entangled state of two atoms in an optical cavity. 

Cavity driving: A scheme for engineered decay mediated by \Ti) 

In the following we briefly discuss a possible |Ti) scheme that combines elements of the 
|To) and of the l^i) scheme. It exhibits an improved error and spectral gap compared 
to the |To) scheme and is also suitable for cavity driving = 0, but not for transversal 
laser driving without transversal confinement. 

Setting 6 = ^ and accordingly, A = g\J^, shifts one of the dressed states of |Ti) and 
|11)|1) into resonance. Then, \S) is effectively prepared through |Ti) by spontaneous 
emission. Similar to the |S'i)-scheme, a choice of /3 = ^^^p guarantees an equal share of 
|T) in the dressed triplet states so that these states decay equally rapidly into \S). 
As compared to the \Tq) scheme the contrast between the unwanted |S'i)-mediated terms 
and the desired |Ti) -mediated terms is more pronounced than previously for the decay 
through 1 5*1) and \Tq). The error and spectral gap are therefore improved compared to 
the I To) scheme: 

(1 - Fs)^T^) = ^ (72) 

The performance of this |Ti) scheme at increased optical driving will also be addressed 
numerically in the section below. 
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7. Comparison of the presented schemes 



In the following we provide an overview of the presented schemes and compare their 
error in the preparation of the maximally-entangled singlet state and their spectral gap 
as a measure for the rate of convergence. We separately discuss the scaling of the static 
error due to the imperfections of the cavity (as discussed in Sec. |4]), and the dynamic 
error originating from increased optical driving by dressing of the levels (as in Sec. [s]). 
An overview of all the schemes and in which section they can be found is shown in Table 
[T] along with a few key results on the performance of each scheme. 



a) 10 



0.5 




10 20 30 40 50 60 70 80 90 100 
Cooperativity C 



10 20 30 40 50 60 70 80 90 100 
Cooperativity C 




0.0 



2.0x10' 4.0x10' 6.0x10' 8.0x10' 
Spectral gap (g) 



Figure 11. Comparison of static and dynamic error for the presented schemes, 



obtained numerically from the full Liouvillian of Eq. (10). (a) Scaling of the fidelity 



with the cooperativity (inset: error weighted with the cooperativity). The lowest error 



is found for the \ Si) scheme (2I;, red solid), followed by the j^o) scheme blue dash- 
dot) and the schemes suitable for cavity driving, |Ti) „„a \^_\ 1 11 



2C ' 

, green dash) and \Tq) (2^, 



sky-blue short dash). Noticeable is the qualitative difference between linear scaling 
of these schemes and the square-root scaling law for the WS scheme { ^J^ ^ black 
dash-dot-dot). Analytical results for the asymptotic scaling are shown for the 
and WS scheme (red and black dots) (b) Fidelity vs. speed of convergence (spectral 
gap). A compromise between fidelity and speed limits the performance of the WS 
scheme, while close-to-linear scaling of both the 151/0) ^-nd |Ti/o) schemes allows rapid 
state preparation. For b) the cavity parameters (7, k) = (^, If ) [32], with C ~ 17 



8 ' 32^ 

were used, in a) C was varied keeping the ratio 7/K = 12/5 constant. The optimized 
parameters used for each of the presented schemes are specified in the corresponding 
section. The same line format is used to denote the schemes in (a) and (b). 



Static error scaling with the cavity parameters 



In Fig. 11 a) we have plotted the fidelity of the steady state with the maximally- 
entangled singlet state, as a function of the cooperativity C = ^ for all schemes 
presented in this work. The curves were obtained numerically by extracting the steady 
state from the full Liouvillian of Eq. (10), using the optimized parameters specified in 
the corresponding section. For the error scaling of the IS*!) and the WS scheme, we plot 
the analytical curves along with the numerical ones. The linear scaling of the static error 
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with C is more clearly seen in the inset, where we plot the weighted error (1 — Fs) C. 
In agreement with our analytic results, we find this quantity to be independent of C 
for C 10 for the |To), |Ti), \Sq) and IS*!) schemes, while the adapted Wang-Schirmer 
(WS) scheme of App. [o] exhibits an error scaling oc \/C as is the case for coherent 
unitary protocol^ The best error scahng of ^ is provided by the IS*!) scheme; an 
increase of the error to ^ is found for the l^o) scheme. The schemes that are suitable 
for cavity driving in the absence of transversal confinement, |Ti) and \Tq) also exhibit 
linear scaling with the cooperativity with further increasing proportionality factors ^ 
and Above cooperativities of C ~ 10 the square-root scaling error of the WS scheme 
is outperformed by the |Ti) scheme which uses similar conditions. An overview 
and numerical examples are given in Table [T] 



Dynamic error scaling with the speed of convergence 

In addition to the static error scaling of Fig. fllla) we present the dynamic error scaling 



with the spectral gap in Fig. 11 b). These curves were obtained by numerically extract- 



ing the spectral gap from the full Liouvillian of Eq. (10). 

Again, the best performance is shown by the l^i) scheme, followed by the other three 
schemes which all exhibit an almost linear scaling; the schemes suitable for cavity driv- 
ing, |Ti) and |To), have a steeper slope. On the other hand, the performance of the 
adapted WS scheme is governed by a compromise between fidelity and speed, that also 
affects the preparation speed. Here, the fidelity with the entangled state drops consid- 
erably at increased speed, so that the state preparation of the WS scheme is found to 
be slow (cf. Ref. [20] )■ Numerical examples of the performance are also given in Table [ll 



We conclude that all relevant benchmarks, both static and dynamic are best for the 
\Si) scheme that was discussed in detail in Sec. [i] and [s] Theoretically, this scheme 
allows for the generation of a maximally-entangled state with fidelities of more than 
90% and convergence time of about 10 /is for present-day cavity experiments. 
Yet, for a possible experimental realization of steady-state entanglement in optical cav- 
ities in the absence of transversal confinement of the atoms, the cavity-driven schemes 
|Ti) and |To) and the randomly laser-driven \Tq) /\Sq) scheme are more suitable. Despite 
the lower proportionality factors, these schemes provide preparation of an entangled 
state with convergence times of a few tens of microseconds at fidelities ~ 80% with 
present day optical cavities. 



8. Conclusion and outlook 



In this article, we have performed a detailed study of the dissipative preparation of 
a highly entangled steady state of two A-atoms in a single-mode optical cavity by 
engineering the naturally occurring sources of noise: spontaneous emission and cavity 

§ See |31j for a detailed discussfon of this point. 
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Table 1. Comparison of the discussed schemes: Analytic scahng of the error and 
rate of convergence (spectral gap) with respect to the cavity parameters for C ^ 10; 
comparative numbers for a cavity QED system as in Ref. |32] , (g, 7, n) /2tt = (16,6, 2.5) 
MHz, C ~ 17. Characteristic dynamic measures are given for a driving that causes a 
dynamic error of 2%. 



Scheme 


Static error 


Spectral gap 


Convergence time 


Transversal confi- 




max. fidelity 


at 2% error 


at 2% error 


nement required? 


\Si) (Sec. g 


3/2 C-i 


172/127 








92.5% 


6 • 10-^5 


10 us 
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\So) (Sec. |ll[ 


7/2 


(5-^5)172/167 
3-10-35 








84.2% 


20 ^ls 
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\Ti) (Sec. § 


9/2 
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60 


no 


\n) (Sec. g 


11/2 


(2 - 73)172/87 
8 • 10-45 








77.2% 


80 /is 


no 


|To)/|5o) (Sec. g 


9/2 


I72/IO7 
1 • 10-^5 








79.7% 


60 /is 


no 


WS (Sec. 


3/2 V2C^^ 


252f72/3A2K 








77.3% 


9 • lO-iy 


70 /IS 


no 



loss. We have employed an effective operator formalism to identify and understand 
the effective decay processes. The schemes we have proposed and analyzed use various 
engineered effective decay processes of either spontaneous emission or cavity loss to 
rapidly reach a maximally entangled singlet state as the steady state of the dissipative 
time evolution at high fidelity. 

Our schemes are suitable for various experimental situations and require coherent driving 
by only a single laser field and another microwave or Raman field; in particular we have 
proposed schemes which work in the absence of trapping of the atoms in the cavity in 
the transversal direction, some of which are tailored for cavity driving. 
We have shown that all our schemes provide a favorable scaling of the static error that 
is linear with the cooperativity of the cavity. In addition we derived the scaling of the 
dynamic error and resolved their underlying mechanisms. Building upon our results 
we have investigated the optimal conditions for the preparation of an entangled steady 
state for a given preparation time. 

We consider our study relevant for the demonstration of an entangled steady state 
by means of dissipation in today's cavity QED experiments. A thorough theoretical 
understanding of the mechanisms allowing for dissipative state preparation of two qubits 
is important as a stepping stone for more complicated studies involving many particles. 
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9. Appendix A: Wang-Schirmer (WS) scheme generalized to A-atoms 

In the foUowing we generahze the scheme of Wang and Schirmer [20J, originaUy proposed 
for two two-level atoms, to A-type atoms. In contrast to the schemes presented so far, 
the ground state |1) of the two atoms is shifted asymmetrically for the two atoms so 
that a coherent coupling is created between IS") and |T). This means that a pure singlet 
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state \S) can, even in the absence of spontaneous emission, no longer be reached as the 
steady state of the time evolution. However, an engineered cavity decay process between 
the triplet states |00) — )■ \T) — )■ |11) is used to prepare a steady state which has a minor 
overlap with |11) and a high fidelity with the singlet. 

By elimination of the excited states we will reduce the coupled A-atom systems to an 
effective system of two coupled qubits described by the master equation of Ref. [20] and, 
subsequently, derive the error scaling of the preparation of the entangled state with the 
cavity parameters. 

The coherent interactions are given by the Hamiltonian of the system as in Eqs. ([6]j9]). 
Here, we use a ground-state Hamiltonian 

4 = E (/3 + (-iy&)|l).(l|+fiMw(|0),(l| + |l),(0|), (74) 

J=l,2 

where a static magnetic field h has been introduced that results in a shift of ground state 
|1) with opposite signs for the two atoms. The driving fields for both atoms exhibit the 
same phase (0 = 0) so that the general effective Lindblad operators are the same as for 
the |T) schemes. The effective Hamiltonian consists of shifts of the ground states 
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1 00) (00 1 - Re 



\S){S\+H,. 



4A 



2, cff 



\T){T\ 



(75) 



A parameter choice of (5 = 0, A ^ (7 ^ (7, k, f2, f^MW; /3, &) and An ^ g"^, allows for 
an adiabatic elimination of the excited atomic levels |e). In this limit, all propagators 
A~gg are simply determined by the shifts of the atomic excited levels A, 



A~V 

We then obtain the effective Lindblad operators 
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(77) 
(78) 
(79) 



where Keff 



A2k 



and 7eff ■■ 



16A2 



The effective Hamiltonian is given by 



eff 



|00)(00| +2/3|ll)(ll| + 



2A 

4A 

b{\S){T\- 



+ /3 
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\T){S\) + ^Mw (|00)(T| + |T)(11| + H.c.) 



(80) 



The corresponding effective couplings are shown in Fig. 12 
In order to match the master equation of Ref. 



we compensate the shifts in the 
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effective Hamiltonian in Eq. (80); i.e. we set the (symmetric) detuning of |1) to 



/3 = — 1^. We then obtain the effective Hamiltonian and the effective cavity decay 



H. 



eff 

TK 



-h {\S){T\ + \T){S\) + r^MW (|00)(T| + |T)(11| + H.c.) (81) 
/^(|ll)(T| + |r)(00|). (82) 

From here we construct the Dicke-type master equation along the lines of Ref. [20]. In 
contrast to Ref. [2DJ, however, we analytically take the (for this scheme) undesired but 
unavoidable spontaneous emission into account and write 



P 



eff, P 



+ P[4\] + ^P[L 



7,fc-| 
eff J 



(83) 



with as the engineered decay and Li as the undesirable spontaneous emission 
processes. In order to analyze the scheme we note that 




b) 




Figure 12. Effective processes of the Wang-Schirmer scheme adapted to A-atoms. 
(a) In the shuffling picture, the triplet states are coupled by a microwave or Raman 
transition fluw and decay by an effective cavity decay KqS- The singlet \S) is 
coherently coupled to \T) by the level shift b. (b) In the steady-state picture, the 
desired steady state 1-05) is no longer coherently coupled, but dissipatively prepared 
from |r) at a rate of KcS,s- Population in li^s) is lost by spontaneous emission j^s- 



1^^^ = /n2 (^1^^) + ^Mwl^)) (84) 

is a steady state of the Hamiltonian of Eq. (77) and also of the Liouvillian of Eq. (83) 
in the absence of spontaneous emission. To understand the dissipative state preparation 
mechanism we change into a basis consisting of {ips) and the orthogonal state 

1^1^ = /n2 ^ u2 (^Mwlll) - b\S)) . (85) 

As can be seen from Fig. 12 b), the singlet-like steady state \ips) is prepared at a rate 
f^e{i,s = \{i's\L'^s\T)\'^ = ^2^^2fe2+r2g,^) ' ^'^''^'^ ^ '^''^'^ decays only by spontaneous 
emission. 

For the derivation of the error scaling we use a rate argument to compare the decay 
rates into and out of the steady state (P = 0, P^^ ~ 1) 

n p ^~QP ^ 3r^g^ 37ft: (4&^ + 3r^Mw) ^Mw /of^N 

(1 - i-^j ~ 3Pt ~ ^f^— p? — ^7:;27^^wV7^^^7;^' ^^^) 
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where we have used the strong couphng condition, and the assumption that the 
populations of the three undesired states are well-shuffled by ^mw so that they have a 
similar population. In contrast to the previously presented schemes, the static error of 
the protocol incorporates a second term that determines the preset compromise in the 
fidelity of the steady state due to the asymmetry h so that 

{l-Fs) ={l-F^^) + {l-\{ijs\S)\') = 

Q4gH2b' + nl^)b' ^2b' + nl^- ^ ' 

The minimal overall error is reached for a trade-off at which these terms are equal. This 
compromise between establishing the steady state by the asymmetry, and at the same 
time avoiding the decrease in its fidelity by the asymmetry, is the cause of the different 
scaling of the error and speed discussed in Sec. [7} For the parameter b we obtain the 
condition 

Inserting this yields the effective decay rate into {ips) 

Ag'n' 2g'n' , , 

«eff,5 = 7 ~ —r^ 89 

^^«(2+7fe) '''' 

and the error of the protocol 

37K (8^ + v^v^) 
4^(37^ + 472^^7^) ' 2V2C' 

Thus, we find that for the WS protocol the fidelity of the steady state with the 
maximally entangled singlet state exhibits a scaling with one over the square root of 
the cooperativity of the cavity. As with the IS*!) scheme, we have one prominent decay 
process to prepare the singlet from the three dressed ground states. With b^pt the 
spectral gap is then given by 

A numerical comparison with the other schemes is given in Sec. [7} 

10. Appendix B: Effects from imperfect coupling of the atoms to the cavity 

Experimental realization of a scheme for dissipative state preparation requires an 
understanding of the effects originating from the imperfect couplings of the atoms to 
the cavity mode. In state-of-the-art cavity QED systems, such as Ref. [32], longitudinal 
confinement prevents fiuctuation of the atomic positions along the cavity axis. Still, a 
static difference in the couplings of the two atoms is possible. Expressing these couplings 
as gi = g{l + a) and g2 = g{l — a) the above analyses can still be carried out for the 
mean coupling of g = \ {gi + g2) ■ From the asymmetry a an additional source of error 
emerges. Below, we exemplarily derive this asymmetry error for the \Si) scheme. In 
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Figure 13. Fidelity under asymmetric coupling to the cavity. Numerical results from 
the full Liouvillian (solid green) are well-approximated by the analytic findings (blue 
dash). High fidelities of about 90% are achieved up to |a| « 0.1, where ^ w 1.22. 



case of a static, asymmetric coupling of the two atoms to the cavity the atom-cavity 
couphng can be written as 

= a{gi\e)i{l\ + g2\e)2{l\) + H.c. 

= ag ((1 + a) |e)i(l| + (1 - a)|e)2(l|) + H.c, (92) 

The asymmetry error affects both the dynamics of the populations and the coherences 
so that we use the effective Liouvillian Ccs to derive the steady state, after having 
excluded other sources of error beforehand (/teg — )■ 0). For weak driving and strong 
coupling g ^ (7, k) ^ {Q, Quw, P) the asymmetry error a can be effectively decoupled 
from both the static error (oc C~^) and the dynamic error (oc VL^) and we obtain 

(1 - Fs)^ ^ 3a'. (93) 



The result is plotted in Fig. 13 using the parameters of Ref. ^2]. The effect of an 
asymmetric coupling is found to be rather small as compared to other sources of errors. 
For |a| ~ 0.1 the loss of fidelity through asymmetry is as little as ~ 2%. In this case, 
with — ^ 1.22, fidelities of about 90% are still achievable. 

11. Appendix C: A scheme for directed cavity decay via l^o) 

In Ref. [22j we presented a scheme for the dissipative preparation of entanglement 
that employs strongly engineered cavity decay to prepare the singlet state IS"). In the 
following we corroborate the claims made about error and speed in Ref. [22]. For this 
1 5*0) scheme, we will drive the atoms with opposite phase = vr so that we can conduct 



our discussion based on the effective operators of Eqs. (31 33) previously derived for 
the \Si) scheme. 



In brief, the mechanism, as visualized in Fig. 14 is the following: Population from 
state 1 00) is driven up to the excited state 15*0) = | (|0e) — |eO)), at a laser detuning 
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Figure 14. Mechanism and effective processes of the j^o) scheme, (a) Setting ^ — ^ 
shifts the lower dressed state of j^o) and |5')|1) into resonance with the optical driving 
17, while the dressed states of |Ti) and |11)|1) are detuned due to a different coupling 
strength of ^/2g. Hence, effective cavity decay k+ from |00) into \S) is enhanced while 
cavity loss k_ from \ S) into is suppressed. The populations of the triplet states are 
shuffled by a microwave or Raman transition IImWi shown together with the effective 
decay processes in (b). 



of A with /3 = 0. State l^o) is in turn coupled by the atom-cavity interaction to 
15)11) with a strength of g. \S)\1) decays into \S) via cavity decay at a rate of k. 

2 

Setting the cavity detuning to 5 = ^ greatly enhances the effective cavity decay 

|00) ^ \Sq) a 15)11) A 15). As in the |To,i) schemes, this is due to the fact that 
the lower dressed state of |5o) and |5)|1) is shifted into resonance. Loss of population 
from the singlet via |Ti) is once again effectively suppressed, since |Ti) and |11)|1) 
are coupled with a larger strength a/25', shifting the dressed states out of resonance. 
The triplet states are shuffled by a microwave or Raman field with optimal strength 

We find that the optimal atomic detuning is A = g\f].- At this detuning we obtain the 
effective operators 



Ll^ =0^|ll)(5|+Zv^|5)(OO| (94) 
LJ^'^^'^^ = ±.v^|00)(00| -^-{\T){S\ ± |5)(5|) (95) 
^S'^'-'^ = -v^|ll)(5| +.V^(±|T)(00| - |5)(00|) , (96) 

where we have assigned «;+ = ^ = 87+, k_ = |^ = 167^. Indeed, the most pronounced 
process is the strongly enhanced effective cavity decay from |00) into \S). The static 
error and the spectral gap are found to be 

(l-^^)|5o)=^ (97) 

\s,) - — - (98) 

A comparative numerical study of the performance of this scheme is given in Sec. [7j 



